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Abstract 

We consider symmetric powers of a graph. In particular, we show that 
the spectra of the symmetric square of strongly regular graphs with the 
same parameters are equal. We also provide some bounds on the spectra 
of the symmetric squares of more general graphs. The connection with 
generic exchange Hamiltonians in quantum mechanics is discussed in an 
appendix. 

1 Introduction 

The symmetric k-th power X^'^^ of a graph X is constructed as follows: its 
vertices are the fc-subsets of V{X), and two fc-subsets are adjacent if and only 
if their symmetric difference is an edge. As an example, and a test case, the 
symmetric square of the complete graph Kn is its line graph. (Useful procedures 
for constructing symmetric squares of arbitrary graphs will be given in Theorem 
4.1 and Lemma 9.1). Symmetric powers were introduced in [1]. 

The symmetric powers are related to a class of random walks, where one 
starts with k particles occupying k distinct vertices of X, and, at each step of 
the walk, a single particle moves to an unoccupied adjacent site. More formally, 
we can generalise the concept of a walk on a graph to a /c-walk, which is an alter- 
nating sequence of /c-subsets of vertices Vi and arcs e^, (Vb, ei, Vi, 62 . . . , e„, Vn), 
such that the symmetric difference of I^-i and Vi is the arc e^. It is readily seen 
that a A:- walk on X corresponds to an ordinary 1-walk on X^''\ 

Our motivation for studying symmetric powers arises from its relevance for 
physically realisable systems and for the graph isomorphism problem. A brief 
outline of the connection between symmetric powers and exchange Hamiltonians 
in quantum mechanics is given in the appendix. 

The relevance to the graph isomorphism problem arises because invariants of 
the symmetric powers of X are invariants of X. There are examples of cospectral 
graphs X and Y such that X^^^ and yf^l are not cospectral. In fact we have 
verified computationally that graphs on at most 10 vertices determined by the 
spectra of their symmetric squares. On the other hand, the main result of this 
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paper is a proof that if X and Y arc cospcctral strongly-regular graphs then 
X^"^^ and y't-^} are cospectral. There is also a family of five regular graphs on 24 
vertices whose symmetric squares are cospectral. Nevertheless, in each of those 
cases, and, in fact, for all graphs wc have examined (including strongly regular 
graphs on up to 36 vertices), the spectrum of the symmetric cubes determine 
the original graphs. (The computations on the the strongly regular graphs on 35 
and 36 vertices were performed by Dumas, Pernet and Saunders; more details 
are given in Section 10.) 

If it were true for some fixed k that any two graphs X and Y are isomorphic 
if and only if their fc-th symmetric powers arc cospcctral, then we would have 
a polynomial-time algorithm for solving the graph isomorphism problem. For 
a pessimist this suggests that, for each fixed k, there should be infinitely many 
pairs of non-isomorphic graphs X and Y such that X^''^ and F't'^J' arc cospectral. 

In the last section of the paper we will consider bounds, from an algebraic 
perspective, on the spectra of the symmetric squares of arbitrary graphs. 

While the focus of this paper is on the spectra of the symmetric squares, it 
should be noted that multivalued graph invariants based on generic (analytic) 
matrix valued functions /(A^'^^) can also be considered [1], where A^''^ is the 
adjacency matrix of X^'^^ In [2] this approach was followed, and numerical 
computations showed that the values of exp{iA^'^^) sufficed to distinguish all 
strongly regular graphs up to around 30 vertices. 

2 Preliminaries 

If A is square matrix, then let 4>{A,t) denote the characteristic polynomial 
det(<I — A) of A. If A is the adjacency matrix of X, we will also write 4){X,t). 
If X and y are vertices of X, we write a; ~ y to denote that x is adjacent to v. 

A graph is strongly regular with parameters (v, k; a, c) if it is not complete 
or empty, has v vertices, and the number of common neighbours of two vertices 
X and y is k, a OT c according as x and y are equal, adjacent, or distinct and 
not adjacent. Thus if X is strongly regular, the neighbourhood of each vertex 
in X is regular and the neighbourhood of each vertex in the complement of X 
is regular. The line graph of the complete graph if„ is strongly regular if n > 3. 

The main tool in this paper will be walk-generating functions. If A is the 
adjacency matrix of the graph X, then the walk-generating function W{X,t) is 
the formal power series 

^A'-f. 

r>0 

We view this either as a power series with coefficients from the ring of matrices, 
or as a matrix whose entries are power series over M. Its zj-entry Wij {X, x) is 
the generating function for the walks in X that start at the vertex i and finish 
at j. 

li D C V{G), then Wd,d{X, t) denotes the submatrix of W{X, t) with rows 
and columns indexed by the vertices in D. The following identities are proved 
in Chapter 4 of [3]. 
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2.1 Theorem. If D is a subset of d vertices of X, then 



t-Uet{WDM^,t~^)) = 



cl>{X,t) 



□ 



2.2 Corollary. Ifi G V{X), then 



t~^Wi,i{X,t-^) = 



(j>{x,t) ■ 



□ 



2.3 Corollary. If i and j are distinct vertices of X 



t-^Wij{X,t-^) = 



{ct>{X\i,t)^{X\ 



j,t)-ct>{X,t)cf>{X\ij,t)y/^ 



□ 



cf>{X,t) 



The presence of the square root in the previous identity is surprising. Note 
though that it causes no ambiguity, since we know that the coefficients of 
Wij{G,t) are non- negative. 

We apply these identities to obtain information about strongly regular graphs. 
If X is strongly regular with parameters {v, k; a, c) and adjacency matrix A then 



(This is essentially the definition of "strongly regular" expressed in linear al- 
gebra.) Since X is regular A and J commute, whence we sec that for each 
non- negative integer k, the power A'^ is a linear combination of I, J and A. 
Thus the generating function Wij {X, t) depends only on whether the vertices 
i and j are equal, adjacent, or distinct and not adjacent. Using the corollaries 
above, this leads to the following: 

2.4 Theorem. Let X be a strongly regular graph. Then (j){X\i, t) is indepen- 
dent of i and, if i ^ j, then (p{X \ ij, t) only depends on whether i and j are 
adjacent or not. □ 

2.5 Theorem. Let X be a strongly regular graph and let Di and be induced 

subgraphs ofV{X). If Di and D2 are cospectral with cospcctnd complements, 
then X\Di and X\D2 are cospectral with cospectral complements. 

Proof Suppose D C V{X). Then WD,DiX,t) is the submatrix oiW{X,t) with 
rows and columns indexed by the vertices in D. Since X is strongly regular, we 
have 



where a, (3 and 7 are generating functions and A{D) is the adjacency matrix of 
the subgraph induced by D. So 



det(al + /?JI + lA{D)) = det((al + 7^(£'))(I + [al + A{D)Y^ pj)) 
= det(al + 7A(D)) det(I + {al + 7A(D))-i/3JJ) 



A^- 



[a - c)A -{k- c)I = cj. 



Wd,b(X, t) = al + /3JJ + 7A(D) 
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Recall that if the matrix products BC and CB are defined then 
det(I + SC) = det(I + CB). 

Since 

= det(I + (al + cA{D))-^hl\^)]> = 11^ 

it follows that 

det(I + (al + -fA{D))-'^pS) = 1 + /3l'^(al + jAiD))-^!. 
We are working effectively over the field of real rational functions in t, therefore 

det(al + jA{D)) = 7!^! det (^^I + A{D)j 

and 

l{aI + jA{D))-'l^ = a-'J2[iy {1,A^1). 

r>0 

We conclude that detWD,D{X,t) is determined by 

a, /?, 7, <l>{D,t) 

and the series 

^{i,A{Dyi)e 

r>0 

which is the generating function for all walks in D. By Exercise 10 in Chapter 
4 of [3] , this generating function is determined by the characteristic polynomial 
of D and its complement. 

Consequently we have shown that if Di and D2 are induced subgraphs of X, 
cospectral with cospectral complements, then X\Di and X\D2 are cospectral. 
Applying this to the complement of X, which is also strongly regular, we deduce 
that the complements of X \ Di and X \ D2 are cospectral. □ 

If 5*1 and S2 are independent sets of the same size in the strongly regular 
graph X, the previous theorem implies that X\Si and X\S2 are cospectral. 
Even this special case of the theorem appears to be new. 

3 Equitable Partitions 

We will also be working with equitable partitions of graphs. A partition n of the 
vertices of X is equitable if for each pair of cells Ci and Cj of tt there is constant 
bij such that each vertex in C'i has exactly bij neighbours in Cj. The quotient 
graph X/tt has the cells of w as its vertices, with bij directed edges from Cj 
to Cj. If G is a group of automorphisms of X, then the orbits of G form an 
equitable partition. If X is strongly regular and u G V{X), the partition with 
three cells consisting of {u}, the neighbours of u, and the vertices at distance 
two from u is equitable. 
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If TT is a partition, the characteristic matrix of tt is the matrix with the 
characteristic vectors of the ceUs of tt as its columns. (Thus it is a 01-matrix 
and each row-sum is equal to 1.) If tt is an equitable partition of X with 
characteristic matrix R and B := A{X/-k), then 

AR = RB. 

There is a matrix B such that AR = RB if and only if col(J?) is ^-invariant, 

and this in turn holds if and only if tt is equitable. If z is an eigenvector for 
B with eigenvalue A, then Rz is an eigenvector for A with eigenvalue A. This 
shows that each eigenvalue of B is an eigenvalue of A. 

As a particularly relevant example, the symmetric square X^"^^ has two sorts 
of vertices: the pairs uv where u v and the pairs uv where u ^ v. If X is 
strongly regular with parameters (w, fc; a, c), then this partition is equitable with 
quotient matrix 

/2a 2A; - 2a - 2\ 
\2c 2fc - 2c )' 

If 5 := a — c, then the eigenvalues of this matrix are 

A; + (5 ± sj{k - (5)2 - 4c, 

and these are eigenvalues of the symmetric square. The eigenvector z oi B 
corresponding to the positive eigenvalue if positive, and therefore Rz is a positive 
eigenvector of A. This implies that the positive eigenvalue is the spectral radius 
of the symmetric square. 

We have the following relation between walks in X and X/tt when tt is 
equitable. 

3.1 Lemma. Let X be a graph with adjacency matrix A. If n is an equitable 
partition of X and B := A{X/tt), then the r,s-entry of B'^ is equal to the 
number of walks of length in X that start on a given vertex in cell Cr and finish 
on a vertex on Cg ■ 

Proof. Assume v = \V{X)\. Let tt be an equitable partition of X with r cells and 
let R be the characteristic matrix of tt. Then AR = RB and, more generally, 

A''R = RB'', k>0. 

Let ei, . . . , e„ denote the standard basis of and let /i, . . . , denote the 
standard basis of W. Let u and v be vertices of X that form singleton cells of 
TT, and suppose {v} is the j-th cell of tt. If u € V{X) then 

{eu,A''Rfj) 

is the number of walks of length A: in X that start at u and finish on a vertex 
in the j-th cell of tt. On the other hand, if vertex u is in the i-th cell of tt, then 
Rcu = fi and 

{eu,RB''fj) = {fi,B''fj). □ 
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4 Constructing the Symmetric Square 



The main result of this paper depends on the observation that we can construct 

the symmetric square of X in two stages. 

We begin with the Cartesian product of X with itself, which has adjacency 
matrix 

A®l + l® A. 

The vertex set of the Cartesian product X nY oi X and Y is V{X) x V{Y), 
and {x, y) ~ (a;', y') if either x = x' and y ^ y' , or x ~ x' and y — y' . We also 
have 

distxny((2:,2;), {x' ,y')) = distxix,x') + disty (y, y'). 
We denote X nX by X^"^. The subgraph of X°^ induced by the vertices 

{{i,i):i€V{X)} 

is called the diagonal. 
The map 

is an automorphism of It fixes each vertex in the diagonal and partitions 

the remaining vertices into pairs. We will call it the Hip automorphism of 

4.1 Theorem. Let X be a graph, let D denote the diagonal of X°'^ and let n 
be the partition of {X^'^)\D formed by the non-trivial orbits of the flip. Then 
is isomorphic to \ D) /tt. □ 

We make some comments on the quotienting involved. Suppose i and j are 

distinct vertices in X. Then ^ (j, and therefore each orbit of the flip 
of size two is an independent set. U i ^ £ and ~ {i,£), then ^ 
Hence two orbits of the flip are either not joined by any edges, or else each 
vertex in one orbit has exactly one orbit in the second. It follows from this that 
((X°^) \ D)/tt has no loops and no multiple edges — it is a simple graph. 

Our aim now is to show that if X and Y are strongly regular graphs with 
the same parameters, then the graphs obtained by deleting the diagonal from 
X°^ and Y nY are cospectral (with cospectral complements). We will then 
show that the quotients modulo the flip are cospectral. 

5 Deleting the Diagonal 

If 6 is an eigenvalue of A, let Eg denote the orthogonal projection onto the 
eigenspace belonging to 6. Then if r > 0, we have the spectral decomposition: 



from which we have 

W{X,t) = Y,{^-iS)~^Ee. 
e 
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Since 

A®I + I(g)A = ^(6* + t)E0 (g) Er, 
e,T 

we see that 

W{X°^, t) = ^(1 - t{0 + r))-'E0 (g) Er. 
e,T 

If M and A'' are mxn matrices, their Schur product (also called Hadamard 
product) M o N is the mxn matrix given by 

{MoN)ij=Mi^jN^j. 

5.1 Theorem. If D denotes the diagonal of X^'^ and A{X) has the spectral 
decomposition J2e ^^o, then 

Proof. It is enough to note that 

iEe(E)E^)D,D= EgoEr. □ 

The linear span of the principal idempotents of the adjacency matrix of a 
strongly regular graph is equal to the span of I, A{X) and J, and is therefore 
closed under the Schur product. Hence Eg o Er is a linear combination of 
principal idempotents. The coefficients in this linear expansion are known as 
the Krein parameters of the strongly regular graph, and are determined by 
the parameters of the graph. Therefore the eigenvalues of Wd,d{X^'^ ,t) are 
determined by the parameters of X, and so det{WD,D{X^'^ ,t)) is determined 
by the parameters of X. 

5.2 Lemma. If X is a strongly regular graph and D is the diagonal of X'-'^, 
then the spectrum of X°'^\D is determined by the spectrum of X. □ 

6 Flipping Quotients 

We use Y to denote the quotient of X°^ by the flip. By Lemma 3.1 we have 
the following. 

6.1 Lemma. IfY denotes the quotient of X'-'^ by the flip and D denotes both 
the diagonal of X'-'^ and the image of D in Y, then 

cl,{Y\D,t) _ <j>{X°\D,t) 

<A(y,t) (l>{x°\t) • 
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Wc now show that, for any graph X, the spectrum of Y is determined by 
the spectrum of X. Given the above lemma it follows immediately that if X is 
strongly regular, then the spectrum of X^^^ is determined by the spectrum of 
X. 

Let Xi and X^ be two cospectral graphs on v vertices with adjacency ma- 
trices A\ and A2. Let L be an orthogonal matrix such that 

L^AiL = A2. 

Let F be the permutation matrix that represents the flip on (gi . So F 
maps x0y to y^x, for all x and y in R". Let R be the normalized characteristic 
matrix of the orbit partition of the flip R is obtained from the characteristic 
matrix of the orbit partition by normalizing each column. We have 

R^R = I, RR^ = ^{I + ¥). 

Let A™ denote the adjacency matrix of X?^. Then there are matrices Ci such 
that 

AfR = RCi 

We prove that Ci and C2 are cospectral. 
We have 

C2 = R^A°^R = R^{L ^ LfAf{L ^ L)R 

whence 

RC2R'^ = RR^{L LfA^^{L L)RR^. 
Because L^L and F commute, so do i (8> i and RR^ . So 

RC2R^ = (L LfRR^A'^^RR^iL i) = (i LfRCi R^{L L) 

and hence 

C2 = R^{L LfRCi R^{L L)R. 

Since 

R^{L LfRR^{L L)R = R^{L Lf{L L)RR^R 

=1, 

we conclude that Ci and C2 are similar matrices. 

Note that it is possible to express the spectrum of Y in terms of the spectrum 
of X. If TT is equitable and B = A{X/tt) and 6 is an eigenvalue of B, then 

dim(ker(B - 91)) = dim(col(i?) n ker{A - 91)) . 

Suppose zi , . . . , z„ is an orthonormal basis for M" consisting of eigenvectors of 
X. Then the products zt zj form an orthonormal basis for R" consisting of 
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eigenvectors of If i ^ j then the span of Zi and zj is equal to the span of 

the symmetric and antisymmetric combinations 

{Zi (g) Zj) + {Zj (g) Zi), {Zi (g) Zj) - {Zj (g Zi) 

These two vectors are orthogonal and the first is constant on the orbit partition 
of the flip, while the second sums to zero on each orbit. If z'^ A = 9z then 
z^RB = 6z^R. So if 6 has multiplicity (. as an eigenvalue of X, the vectors 

[Zi® Zj) + {Zj (g> Zi), Zi®Zi, 

where Zj G ker(A — 91), give rise to a subspace of eigenvectors of Y with eigen- 
value 26 and dimension {^^)- If has multiplicity i and r has multiplicity m, 
then we obtain a subspace of eigenvectors of the quotient with dimension £m. 
By adding up the dimensions of these subspaces, we find that the images of the 
given vectors provide a basis consisting of eigenvectors oiY. It follows that the 
multiplicities of the eigenvalues of Y are determined by the eigenvalues of X 
and their multiplicities. (If X has exactly r distinct eigenvalues, then X^"^ has 
at most ('^2 ^) ; if X°'^ has fewer eigenvalues, then the procedure just described 
will give the multiplicities of the eigenvalues of Y, but does not lead to a simple 
formula.) 

7 More Cospectral 

We have seen that if X and Y are strongly regular graphs with the same pa- 
rameters, then their symmetric squares are cospectral. Here we extend this. 

7.1 Lemma. If X and Y are strongly regular graphs with the same parmeters, 
then the complements of their symmetric squares are cospectral. 

Proof. I Prom Exercise 22 in Chapter 2 of [3], we have 

(l){X,t + l) = - l'^{tI + A)-^l). 

I From this it follows that cospectral graphs X and Y have l cospectral com- 
plements if and only if the generating function for all walks in X is equal to the 
corresponding generating function for Y. 

Assume X is strongly regular, let A denote the adjacency matrix of X^^^, let 

TT be the partition of the vertices of X^"^^ by valency and let the characteristic 
matrix R and quotient matrix B be defined as in Section 3. Then AR = RB 
and so, for if ^ > 0, 

A^R = RB^. 
Since the columns of R sum to 1, 

I'^A^l^ = l^A^Rl = l^RB^l. 

We have 

\^R = {vk/2, v{v-l- k)/2) 
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and therefore the entries of R"'" are determined by £ and the parameters of 
X. Hence the generating function for ah walk in X^"^^ is determined by the 
parameters of the strongly regular graph X, and the result follows. □ 



8 Variations 

The direct product X x Y oi graphs X and Y has vertex set equal to V{X) x 
V{Y), and (u, v) ~ {x, y) if and only if u ~ a; and w ~ y. We have 

A{X X y) = A{X)®A{Y). 

The flip map 

{x,y) ^ {y,x) 

is again an automorphism of X x X that fixes the diagonal. We can obtain an 
analog of the symmetric product by deleting the diagonal and then quotienting 
over the flip. A slightly modified version of the argument in this paper shows 
that if X is strongly regular, then the spectrum of this analog is determined by 
the spectrum of X. The key step is to verify the following analog of Theorem 5.1: 

WD,D(x®^^) = ^(1 - t0T)-i£;e o i;^. 

For a second analog, we turn to the graph obtained from the Cartesian 
power X^'' by deleting the diagonal and the quotienting over the orbits of 
the automorphism that sends each fc-tuple to its right cyclic shift. Again our 
argument shows that if X is strongly regular, the spectrum of this analog is 
determined by X. Thus there is more than one candidate for the "symmetric 
cube" of a graph, but the spectrum of the one just described is a less useful 
graph invariant than the spectrum of the symmetric cube defined in Section 1. 



9 Symmetric Squares of General Graphs 

In this section we take a closer look at the purely algebraic properties of the 
symmetric powers, and of the symmetric square in particular. We start by 
giving a purely algebraic definition. 

Let P(*=) be the 0/1-matrix with Q rows, labelled by the fc-tuples 
with 1 < i < 7 < . . . < ^ < w, and w'^ columns, labelled by the fc-tuples 
[i' ,j' , ■ ■ ■ with 1 < i', j', . . . , r < u, such that the elements P^^\ j^, j, 
are 1 iS ... ,1) is a permutation of I'] . Then 

9.1 Lemma. The adjacency matrix A^''^X) of X^''^ is 

^^'=>(X) = (^PW {A{X)(^lf''-^) P«* 
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Wc focus on the symmetric square, and more generally on the properties of 
the linear map 

: G n{G) = G^'^'i = P(2) (G (8) I)P^^)* • 

Henceforth, we will write P instead of P^^^ . 

Because fl is the composition of the two completely positive maps [6] A t-^ 
A^l and A i— > BAB*, O is completely positive itself. In particular, Cl preserves 
positive semi-definiteness. One easily checks 

PP* = 2I(-d^ (1) 

d V 

P*P = ^{Eii^Ejj+Eij^Eji)-2^Eii^Eii, (2) 

i,j=l i=l 

where {Eij} is the standard matrix basis. 

The spectrum of a general Hermitian matrix and the spectrum of its sym- 
metric square have the same average value. When G is an adjacency matrix this 
obviously has no import, because adjacency matrices are traceless. However, in 
certain quantum mechanical contexts the map Cl is applied to Hamiltonians 
which are not traceless. 

9.2 Theorem. For G av x v Hermitian matrix, 

Tr[G]/t; = Tr[G{^>]/(^ 

Proof. The partial trace of P*P over the second tensor factor, defined as 
Tr[(X I)A] = Tv[X Tr2[A]], yields 

V V 

V V 

= ^ {Eii + EijSij) — 2 ^ Eii 

— l i=l 



= {V-I)lv. 



Therefore, 



Tr[n(G)] = Tr[P*P(G0l)] 
= Tr[G Tra [P*P]] 
= {v-l)Tr[G\. 

Dividing by v{v — 1) yields the statement of the Theorem. □ 
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9.1 Comparison between the spectrum of a matrix and 
the spectrum of its symmetric square 

For a Hcrmitian matrix A, we denote by A^(j4) its k-th largest eigenvalue, 
counting multiplicities. Likewise, A|.(yl) is its k-th smallest eigenvalue. 
We prove the following: 

9.3 Theorem. For any non-negative positive semi-definite v xv matrix G, the 
following relation holds, for 1 <m<v: 

Proof. We refer to [4] or [5] for the basic matrix analytical concepts and 
theorems. 

say it at all :-) Focusing on a particular value of m, 1 < m < we need to 
show 

for all G > 0, or, equivalently, 

AJ„(P(G®I)P*)>1, (3) 

for all G > with X^{G) = 1. 

First note that one needs to prove this only for G a partial isometry of rank 
m. Indeed, for every G > with Xl^{G) = 1, there exists a partial isometry B of 
rank m such that G > B. As noted above, f2 is a completely positive map, hence 
n{G) > n{B). By Wcyl monotonicity we then have X^{Q{G)) > A^(n(B)). 
Thus (3) follows for G if it holds for B. 

Let us write B as B = Q*Q, with Q e Mm^„(C) and QQ* = !„. Let qj be 
the j-th column of Q. Thus the qj are v fc-dimensional vectors and 

V 

The matrix P{Q*Q (g) I)P* has the same non-zero eigenvalues as 

(Q0I)P*P(Q* 01). 
Using the explicit form (2), a short calculation shows that 

xUm*Q)) = XU^ + A) = 1 + Xi,{A), 
where Ais av x v block matrix with blocks Aij of size m x m given by 

Aij = (1 - 26ij)qjq*. 

We have 

i=l 
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Wc have to show that A4i(A) > 0. To that purpose, consider the principal 
submatrix A' of A consisting of the 2x2 upper left blocks: 



A' = 



Au Ai2 
A21 A22 





If we can prove that Xi,{A') > 0, this implies Xi,{A) > via eigenvalue inter- 
lacing. 

When m = 1, the qi arc scalars, and direct calculation shows that A|(^') = 0. 

For m > 1, consider a (non-orthogonal) basis of in which qi and (72 are 
the first basis vectors. Let S be the transformation from this new basis to the 
standard basis. Under the *congruence governed by S, A' is transformed to 

/ 



SA'S* 



-1 







\ 







1 














1 . 


. 











-1 










/ 



V 

This matrix has eigenvalues —1, with multiplicity 3, 0, with multiplicity 2m — 4, 

and 1, with multiplicity 1. By Sylvester's Law of Inertia, a *congruence does not 
change the sign of the eigenvalues. Thus A' has 2m — 3 non-negative eigenvalues 
as well. Hence, for m > 2, X!^{A') > 0. 

To cover the remaining case of m = 2. we first perform a specific * congruence 
on A directly. For m = 2 there are only 2 independent vectors qj. Let 5*1 be 
the transformation that brings qi to (1,0), and 92 to (0,1). Let qs be brought 
to (x, y). We can assume without loss of generality that 53 ^ q2, so that x ^ 0. 
The 3x3 upper left blocks of SiASl will thus be 

/ 














X 


\ 








1 





!J 








1 











X 











-1 





y 


X* 


y* 








-\x\' 


-xy* 








X* 


y* 


-x*y 


-\y\' J 



One further *congruence S2 = l + Ei^5/x* brings this to S2S1ASIS2, with 3x3 
upper left blocks 

/ 






{y/xT 











-x{y/x)* 


y/x 





1 





y 








1 











X 











-1 





y 





y* 










-xy* 


-x*y/x 





x* 


y* 


-x*y 


-\y\' 
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The upper left 3x3 principal submatrix is of the form 




which has eigenvalues and ±-\/l + |zp, i.e. it has two non- negative eigenvalues. 
By the interlacing theorem, S2S1ASIS2 must then also have at least two non- 
negative eigenvalues, and by Sylvester's Law of Inertia, A itself too. □ 

Because of the restriction to positive semi-definite matrices, Tlicorem 9.3 
can only be applied directly to graph invariants formed from, say, the spectrum 
of the Laplacian matrix L{X) of the graph under the map CI. The following 
Corollary extends Theorem 9.3 to Hermitian G that are not necessarily positive 
semi-definite, and can therefore be applied to adjacency matrices proper: 

9.4 Corollary. For any Hermitian v x v matrix G, 

(Ai(G) + Ai(G))/2 < Ai(G{2}/2) (4) 
{Xl{G) + Xl{G))/2 > Al(G{2}/2). (5) 

Proof. Let a = Xi{G), then G' := G + al > 0. Applying Theorem 9.3 to G' 
gives 

XUG + al)<XimG + aI)). 

Noting that f2(I) = PP*, which has the same non-zero eigenvalues as P*P = 
2I^vy yields 

XUG)+a<XU^{G))+2a, 

and the first inequality of the Corollary follows. The second inequality follows 
by applying the first one to — G. □ 

Very likely, the bound of Theorem 9.3 (and the Corollary) can be sharpened. 
However, it cannot be sharpened by more than a factor of 2. This can be seen by 
taking as G a rank- A: partial isometry, for which Xj,(G) = 1, and noting that by 
inequality (5) (with k = v), G < I implies G^^^ < 21. Hence, for this particular 
G, A^(G^^^) < 2A^(G), which would contradict a sharpening of Theorem 9.3 by 
a factor of more than 2. 



9.2 On the nature of P^^^ 

In this section we consider the P*^*) appearing in the definition of the symmetric 
power, and compare it to the two related operators Py and P^, which are 
projections from the fc-fold tensor power of C" to its totally symmetric and 
totally antisymmetric subspace, respectively ([5], Section 1.5). Formally, Pv and 
Pa are defined as those linear operators that map a tensor product of k vectors 
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from to their symmetric and antisymmetric tensor product, respectively, 
Pv{xi ^ ■ ■ ■ Xk) = (g)...(8)a;<^(fc) 

(7 

P/\{xi ■ ■ ■ ^ Xk) = (A:!)"^^e^a;^(i)(8)...(8)a;<^(fe), 

cr 

where the sum is over all permutations a oi k objects, and is the signature 
of a. The operator P^*^^ is similar to Py in that tensor products that differ 
in the ordering of factors only are mapped to one and the same vector; it is 
similar to Pa in that it maps to a space of the same dimension as the totally 
antisymmetric subspace and maps tensor products containing identical factors 
to 0. 

To describe this in a more formal manner, consider the basis of the totally 
antisymmetric subspace consisting of the vectors 

e(i,j,...,o = Cj A A • • • A e; 

:= ^ e^e<,(i) e^^j^ (g) . . . (g) e^^i) , 

a 

labelled by the fc-tuples (z, j ...I) with l<i<j<...<l<d. Then P^''^ 
maps the vector ej' (g) e^' (8) • • • (g) ei>, where [«', j', •••,/'] is a fc-tuple with 1 < 
< d, to the vector e(ij-...^;), with fc-tuplc equal to the 

fc-tuple •■•,?'] sorted in ascending order, provided ■ ■ ■ ,1'] does not 

contain equal indices, and to otherwise. The difference between P^*^' and Pa 
is the absence of the sign eo- of the permutation that realises the sorting. Note, 
for k = 2, 

P;Pa = (I-F)/2 
P*Pv = (I + F)/2, 

where F is the flip operator defined in section 6. 

In the following we look at the map G i-> := Pv(G 1®'=~^)Pv. Because 
of the symmetry of Pv , 

= ^Pv(G(8)l0...(8>I + I(8)G(8)I...(8>I + ...+I(8>I(8)...(8)A)P* 

= H Pv(I + iG)«'^-p*. 
k at t=o 

The expression Pv(I + tG)®'^P* is nothing but the totally symmetric irreducible 
representation of I + tG on fc copies of C. It is well-known from representation 
theory that the eigenvalues of an irreducible representation of a matrix A depend 
only on the eigenvalues of A itself. Therefore, we find that the spectrum of 
P\/{G®l®^~^)Py depends on the spectrum of G only. In other words, if Gi and 
G2 are cospectral, then so are G^ and G2 . A similar reasoning applies when 
using Pa instead of Py. 
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It is therefore remarkable that fl{Gi) and il{G2) need not be cospectral even 
if Gi and G2 are, given that PC^) is a combination of Py and Pa- This is one 
the underlying reasons why we chose to study fl in the context of the graph 
isomorphism, the other reason being its physical relevance (as discussed in the 
appendix) . 

10 Computational Results 

Strongly regular graphs, and to a somewhat lesser extent walk-regular graphs, 
satisfy very strong combinatorial and algebraic regularity conditions, and it 
might be hoped that this was closely related to the occurrence of cospectral 
symmetric squares. Unfortunately our computational results show that this is 
not the case, and that in fact graphs with cospectral symmetric squares occur 
in relative abundance. Nevertheless, the examples that we have found do have 
some interesting algebraic properties that may go some way towards explaining 
when symmetric squares are cospectral. 

We have checked all graphs on up to 10 vertices without finding any pairs 
of graphs with cospectral symmetric squares, and currently the smallest pairs 
that we know have 16 vertices. There are only two pairs of cospectral strongly 
regular graphs on 16 vertices, but using a variety of heuristic search techniques, 
we have constructed more than 30000 further graphs on 16 vertices that have 
a partner with a cospectral symmetric square. These heuristics involve first 
using direct searches of catalogues of strongly regular graphs, vertex-transitive 
graphs and regular graphs to generate an initial collection of example pairs. 
Then we construct large numbers of closely-related graphs by making a variety 
of minor modifications to these initial graphs, such as exchanging pairs of edges, 
removing one or more vertices, removing one or more edges, or adding or deleting 
one-factors. These graphs are then searched for further non-isomorphic pairs 
of graphs with cospectral squares, and any new examples added to the growing 
list. By repeatedly applying these techniques, we can obtain pairs of graphs that 
are seemingly very different to the initial examples, but that have cospectral 
symmetric sqiiarcs. 

Using these techniques, we have found it easy to construct many pairs of 
graphs on 16 or more vertices cospectral squares. We have put considerable 
effort in constructing as many graphs as possible on 16 vertices, but due to the 
techniques involved, we do not speculate as to whether these 30000-1- graphs 
might comprise most of, or almost none of, the full collection of examples on 16 
vertices. All our efforts to construct examples on fewer than 16 vertices have 
failed. 

The examples that we have constructed do not show any strong graph- 
theoretical structure, most of them are not regular, and there are many examples 
with trivial automorphism group. However the pairs of graphs with cospectral 
symmetric squares do exhibit interesting algebraic behaviour that is not a priori 

necessary in order to have cospectral symmetric squares. In particular, for all 
of the known pairs of graphs {X, Y} such that X^"^^ and Y^^^ are cospectral. 



16 



the following properties also hold: 

(a) X and Y are cospectral, and X and Y are cospectral, 

(b) The symmetric squares of X and Y are cospectral, 

(c) The complements of the symmetric squares of X and Y are cospectral 

(d) The multisets {(piX\i) : i e V{X)} and y{Y\i) : i £ V{Y)} are equal, 

(e) The multisets {'p{X\ij) : i,j G V{X)} and y{Y\ij) : ij e V(Y)} are 
equal. 

If X and Y are strongly regular graphs with the same parameters, then all 

of these five properties hold (the third one requires a non-trivial argument) , but 
in general we do not know whether or not these are necessary conditions for X 
and Y to have cospectral symmetric squares. 

There are 32548 strongly regular graphs with parameters (36,15,6,6) each 
of whose symmetric cubes has 7140 vertices. Performing exact calculations of 
characteristic polynomials on matrices of this size requires highly specialized 
software, and the only such software of which we are aware is that being devel- 
oped by the LinBox team (see www.linalg.org). Proving that two graphs are 
not cospectral is easier in that if there is some a G GF{p) (where p is a large 
prime) such that det(Ai + al) ^ det(^2 + ctl) (mod p) then Ai and A2 are 
definitely not cospectral. We would like to thank the LinBox team, particularly 
Jean-Guillaume Dumas, Clement Pernet and David Saunders for planning and 
performing computations using this technique that demonstrated that none of 
the SRGs on 35 or 36 vertices have cospectral symmetric cubes. 
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12 Appendix: Quantum Hamiltonians and Sym- 
metric Powers 

Consider a generic set of n distinguishable two-dimensional quantum systems 
(qubits). Letting |0), |1) be a basis for C^, and defining raising and lowering 
operators for qubit i: 

5+ = |l)(0Ur = |0)(l|, 

a commonly encountered interaction Hamiltonian for the systems is of exchange 
form: 

where gij is the interaction energy between qubits i and j. For instance, the 
systems could be two- level atoms in a molecule, interacting via a dipole-dipole 

interaction; spins on a lattice interacting via an "XK" spin-exchange interac- 
tion; or hard-core bosons hopping around some lattice structure (Bose-Hubbard 
model). 

In certain situations the relevant physics lies only in the properties of this 
interaction Hamiltonian. For instance, for the two-level atoms the free Hamil- 
tonian is trivial and can be ignored by going to the 'interaction picture'. In the 
limit of hard-core bosons in a Hubbard model, the interaction energy dominates 
the single-site energy, and double occupancy of a site is forbidden. In such sce- 
narios, if it is also approximately true that the interaction strength is the same 
regardless of the pair of systems under consideration (no distance dependent in- 
teractions for instance) then we can take g^j = 1,0 according to whether qubits 
i and j are coupled or not. This simplified interaction Hamiltonian is then 

n 
fe=l 

i.e., a direct sum of the symmetric powers of the underlying graph X, whose 
adjacency matrix is gij). 

There arc two main types of graphs that generally come under consideration 
in physics, neither of which are particularly interesting from the graph theoretic 
point of view: (i) Small, (generally planar) graphs corresponding to molecular 
systems. (Does the excitation spectrum of a molecule determine its structure?) 
(ii) Large 'local' graphs in M^'^'^ corresponding to nearest neighbour interac- 
tions - in general some sort of standard lattice structure. In the latter case 
the interesting physical properties (phase transitions, super conductivity, etc.) 
generally appear for a number of excitations k ~ n/2. 
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To understand the strength of graph invariants formed from such Hamil- 
tonians, and the complexity of deahng with such Hamiltonians in physics, the 
following observation (discussed formally in section 9.2) is useful: The subspace 
of the full Hilbcrt space in which the fc'th excitation block of the Hamiltonian 
lives is one of both bosonic and fermionic nature. Although the Hamiltonian 
is strictly speaking bosonic, fermionic features arise due to it not being possi- 
ble for two excitations to reside in the same qubit. Thus, the bosons, instead 
of living in the ("^^~^) dimensional symmetric tensor power subspace V'^H, 
rather live in an "unsigned" version of the antisymmetric tensor power space 
/\^TL. ("Unsigned" refers to the fact that the antisymmetry is not present). 
If, instead of living in such a hybrid "Fermi-Bose" subspace of Hilbert space, 
the excitations were to live in these more standard subspaces, it is easy to see 
that their spectra would essentially be equivalent to that of the single particle 
spectra (the standard graph spectrum). 

Finally, it should be noted that an efficient quantum circuit simulating evo- 
lution under iJjnt is guaranteed to exist by various standard results in the theory 
of quantum computation. This opens up the interesting possibility that graph 
invariants based on symmetric fc-th powers of a graph for k = 0{v) are quantum 
computationally tractable, whereas classical tractability would seem to require 
that k = 0(1). 
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